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We perform a model-independent study of inelastic dark matter at the LHC, concentrating on 
the parameter space with the mass splitting between the excited and ground states of dark matter 
above a few hundred MeV, where the direct detection experiments are unlikely to explore. The 
generic signatures of inelastic dark matter at the LHC are displaced pions together with a mono jet 
plus missing energy, and can be tested at the 7 TeV LHC. 
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Introduction. Despite overwhelming evidence from 
astrophysical observation, we still don't know the parti- 
cle properties of dark matter (DM) or how it interacts 
with standard model (SM) particles. Dark matter as a 
weakly interacting massive particle (WIMP) is predicted 
in many extensions of the SM and motivates searches for 
its direct detection, where one looks for ambient WIMPs 
scattering with heavy nuclei. An implicit assumption in 
these searches is that DM particles can scatter elastically 
off detector nuclei. So far, there are no unambiguous sig- 
nals for a detection of DM scattering, leading to strong 
bounds on WIMPs with weak scale interactions [l[ . 

These null results raise severe questions about the vi- 
ability of the WIMP paradigm, and motivate theoreti- 
cal exploration of modifications of the standard picture. 
For example, if dark matter is required to scatter in- 
elastically into a state heavier than the initial WIMP by 

> 1 MeV, typical WIMPs in the Milky Way halo will 
have insufficient energy to upscatter, explaining the null 
results of direct detection experiments. Nonetheless, we 
shall see below that inelastic dark matter (iDM) mod- 
els can be tested at the Large Hadron Collider (LHC). 
The iDM models were introduced some time ago @, y] 
to evade constraints on models from null DM direct de- 
tection searches, and were recently revisited to reconcile 
the DAMA observation of annual modulation [3| with 
other null direct searches [5j , provided the mass splitting 
is below a few hundred keV. 

In this letter, we explore the LHC's capability to iden- 
tify iDM scenarios through a model-independent, effec- 
tive field theory (EFT) approach [6Hl0j. We discuss 
general strategies and characterize signatures appropri- 
ate for different parameter space. The mass splitting A 
and coefficients of the various EFT operators control the 
lifetime of the excited DM particle, but it is generically 
long-lived on collider time scales. If the mass splitting is 
too small, the decay products of the excited state are too 
soft to be observed and the generic signature is a monojet 
plus missing energy However, for mass splittings 

> I GeV, iDM results in a spectacular signature with 
displaced pions appearing on top of the monojet signa- 
ture. 

Operators and interactions. We assume the dark 
sector is composed of two SM gauge singlet fermions x 



and x* with masses m and m», where A = m» — m > 
characterizes the splitting between the two states. We 
assume that interactions with the SM are required to 
contain one x an d X* each, and (for simplicity) restrict 
ourselves to interactions consisting of operators which 
involve the up-quark and preserve parity: 
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where A, parameterizes the strength of each interaction. 
The effective Lagrangian consists of the SM plus kinetic 
terms for x an d X*i as wen as these operators plus their 
complex conjugates. 

We are interested in mass splittings ranging from 100 
MeV to several GeV. For A < 1 GeV, we can use the 
chiral Lagrangian to describe how these interactions lead 
to x* decaying into x pl us a number of pions. We focus 
on non-vanishing terms with the minimum number of 
pion legs, since higher order interactions are phase-space 
suppressed. The leading terms for each operator are 
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where F n = 184 MeV and (uu) = -(242 MeV) 3 (eval- 
uated at 1 GeV) [H|). For on-shell x an d X*i 01 can 
be rewritten in terms of m and , in analogy with the 
usual treatment of pion-nucleon interactions. The parity 
of the quark bilinears involved dictate that the leading 
terms induced by 0i,2 in the chiral Lagrangian corre- 
spond to one pion interactions, whereas 03,4 correspond 
to interactions involving at least two pions. For theories 
in which dark matter also couples to down-quarks, the 
right-hand side of Eq. @ for O1.4 project onto the iso- 
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FIG. 1: Decay lengths of the excited dark matter state at FIG. 2: Production cross section of pp -> \z Xg j at the 7 TeV 
rest as a function of mass splitting for different operators with LHC. 
Aj = 1 TeV and a dark matter mass of 5 GeV. 



spin violating terms, whereas for 02,3 arise from iso-spin 
conserving pieces. 

Decay of the x*- We concentrate on the regions of 
parameter space with A > m,o (for 0\^) or A > 2m n + 
(for 03,4), for which decays lead to real pions and re- 
sult in hard enough decay products so as to register in 
LHC detectors. We deal with the regimes A < 1 GeV 
(described by the chiral Lagangrian) and A > 1.5 GeV 
(described by interactions with quarks) separately. 

For A < 1 GeV, 0\ and 02, lead to two-body decays 
X* — > X + n ° w ith decay widths r, (i=l,2) 
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where m = (m* + m)/2, a.\ — F^m/Af, and «2 = 
I (tin) I / (F^A 2 ). In the limit m ^> A > m^o, Ti is roughly 
independent of the DM mass, and T 2 oc A 3 /m 2 . Opera- 
tors 03,4 result in three-body decays of x* to x plus two 
pions. Neglecting pion masses and in the limit m 3> A, 



r 3 (x* -> x^+O = 2r 3 ( x * -> x^V) 

A 5 



(uu) 2 A 3 
8^Af ' 



T4(X* -> X7T+7T ) = 



240tt 3 A| ' 



(4) 



For A J> 1.5 GeV, the chiral Lagrangian is no longer a 
suitable description, and we compute x* ~~ ^ X99? 
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where a\ = 1/20, 02 = A 2 /(560 m 2 ), and 03 = 
a 4 = 1/60. The decay produces soft jets of hadrons 
(mostly pions with a small fraction of kaons) described 
by the parton shower of QCD. The intermediate region of 
1 GeV < A < 1.5 GeV is complicated, and receives con- 
tributions from resonances as well as multi-pion states. 



We have approximated the behavior in this region by re- 
quiring the decay length smoothly interpolate between 
the leading decay for small A and the soft jet regime of 
large A. 

Decay lengths cr as a function of A, for Aj = 1 TeV 
and m = 5 GeV are shown in Fig. [TJ Different values 
of Aj rescale the presented lifetime by (A new /1 TeV) 4 . 
For the chosen parameters, they vary from 1 cm to 10 m, 
depending on A and the operator mediating the decay. 
As we shall see below, the most useful production regime 
at a hadron collider results in relativistic x* whose life- 
times in the detector frame are given by cr = jctq, where 
7 = E Xt Im*. In our detailed LHC calculations below, we 
include this dilation factor on an event- by- event basis. 

Production of iDM particles at the LHC. At the 
LHC, the interactions with quarks will result in events 
containing one x an d one x* m the final state. The 
hadrons from the x* decays are generally too soft to 
be used as triggers. However, a "monojet" process, 
PP X* X3 (phis the conjugate XX* j process) contain- 
ing an additional unflavored jet j radiated from the initial 
partons can provide a suitable trigger. After the x* de- 
cay, the final state consists of 2x + j + 7r's, where because 
of the long x* lifetime, the 7r's are produced far from the 
primary interaction vertex, leading to a signature of a 
monojet plus displaced pions. 

Current LHC monojet searches [lU, [l3[ rely on a miss- 
ing energy trigger and apply a missing energy cut of 
# T > 150 GeV (at CMS). For the iDM signature, one 
may use the same triggers and missing energy cut. We 
simulate the expected production rate after the missing 
energy cut at the 7 TeV LHC using Madgraph 5 [3] with 
the CTEQ 6L1 fl5| parton distribution functions (PDFs). 
The results for the various operators are shown in Fig. [5] 
For DM masses below around 50 GeV, the rates become 
independent of the value of the mass itself, because the 
JfiT cut becomes the limiting factor for production. 

In Fig. 02 we show the distribution of the 7 factor for 
the excited WIMP (after cuts) resulting from production 
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FIG. 3: The distribution of 7 (for Oi) for three choices of DM 
masses after a missing energy cut oHfr > 150 GeV. 
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FIG. 4: The transverse energy of hadrons (resulting from op- 
erator 0i ) for two choices of DM masses and splittings, as 
indicated. 



through operator 0\. Results from O2-4 are very sim- 
ilar. As expected, the peak of the distribution shifts to 
larger values for smaller DM masses, and remains ~ 10 
for masses as large as 100 GeV. Comparing Figs. [T] and 
[3] reveals that for WIMP masses around 50 GeV, the lab 
frame decay length jctq is around 1 m for a wide range 
of A and A - 1 TeV. 

Search strategy and discovery potential. In ad- 
dition to the primary jet against which the XX* system 
recoils, iDM also produces one or more hadrons deep in 
the detector from the x* decay. For decay lengths on 
scales of 10 cm to 1 m, the hadrons are likely to deposit 
most of their energy in the electromagnetic (ECAL) or 
hadronic (HCAL) calorimeters. For the discussion be- 
low, we use detailed numbers for the CMS detector [r|, 
though similar conclusions will also hold for ATLAS. The 
hadrons from the x* decay are not typical of QCD jets, 
since they contain a much smaller multiplicity of charged 
particles (particularly when charged pions are produced 
after x* traverses the tracker). Aside from their relatively 
long decay length, these features are common between x* 
and hadronic tau decays (the actual hadron multiplicities 
from x* decays are different from hadronic taus and de- 
pend on the mass splitting). 

For a sufficiently loose definition of a hadronic tau, the 
iDM decays can mimic an extra jet consistent with the 
hadronic tau signature. After being tagged as a hadronic 
tau, one can further reduce real tau backgrounds using 
the decay length. Our search strategy is thus a mono- 
jet together with a tau-tagged jet. At CMS, hadronic 
tau's resulting in visible transverse energy > 15 GeV 
have a tagging efficiency of ~ 20%, with a fake rate of 
0.3% [l7j]- We adopt these numbers as estimates for the 
(mis)tag rates for the x* decayed jet, though it would 
be worthwhile to have a proper treatment by the exper- 
imental collaborations using a realistic detector simula- 
tion. Given the very short r lifetime (cry ~ 100 fim in 
the r rest frame), we estimate that the background from 
Wj, where W — > tv will be essentially entirely removed 



by a displacement cut > 10 cm. 

In the limit of the missing Ifix > ffi and A 3> m n o , the 
total transverse energy of the displaced pions is pt(tt°) ~ 
Et A /m x . Lighter x masses and larger mass splittings 
result in more deposited energy of displaced pions in the 
calorimeter. In Fig. [4] we show pr of ir° (for Oy\ other 
operators lead to similar results) for a few choices of m 
and A. Requiring Pt(^°) > 15 GeV, the signal efficiency 
drops rapidly as the mass varies from m = 50 GeV to 
m = 10 GeV. 

In order to assess the discovery potential at the LHC, 
we allow for events which pass a monojet + "displaced 
hadronic tau" selection. Thus, we select events with: 
Nj = 2 with p T (ji) > 110 GeV and 15 < pr(j 2 ) < 
30 GeV as well as # T > 150 GeV. While in principle 
one could allow the second jet (typically from the x* 
decay) to have a higher px, in practice for these choices 
of A and pt(ji), the signal is not very sensitive to the 
upper bound on prti?)- We require the x* decay occur 
before the barrel ECAL, 129 cm from the center along the 
radial direction [l6[ . We apply a 20% hadronic r-tagging 
efficiency to signal events satisfying these criteria. 

We take the backgrounds measured as part of the stan- 
dard CMS monojet search based on 36 pb _1 [11] and ap- 
ply the 0.3% hadronic tau mistag rate. In practice, this is 
an over-estimate of the background level because the pre- 
sented backgrounds include both one and two jet events; 
however it can only over-estimate the background, and 
suffices for a conservative estimate of the LHC sensitiv- 
ity. To estimate future reach, we rescale the background 
to a 5 fb _1 data set in order to determine the future sta- 
tistical uncertainty (which is still expectd to be the dom- 
inant contribution), thus determining future 95% limits 
on a putative signal at a 7 TeV LHC which has collected 
5 ft.- 1 . 

In Figs. [5] and [5] we show the expected reach in terms 
of the bound on the strength of the contact interaction 
Aj for the axial-vector interaction 0\ and scalar inter- 
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FIG. 5: The 95% C.L. exclusion limit on Ai (for Oi) as a 
function of m, assuming 5 fb _1 of 7 TeV LHC data, for the 
three values of mass splitting A shown. The approximately 
straight lines are limits derived using monojet only searches 
at 36 pb" 1 (solid) and 5 fb _1 (dotted). 
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FIG. 6: The same as Fig. [5] but the scalar interaction O3. 



action O3, as a function of m and three choices of A, 
for the LHC operating at 7 TeV and having collected 5 
fb -1 . Also shown for reference are the limits from stan- 
dard monojet searches from the existing 36 pb" 1 and 
projected with 5 fb -1 . For mass splittings > 1 GeV, 
limits from the displaced pion search provide more strin- 
gent limits on iDM than the standard monojet search. 
For smaller mass splittings, the visible x* decay prod- 
ucts have difficulty passing the px > 15 GeV cut. We 
find that limits on 0\ and O3 are comparable. We ex- 
pect that the vector operator C4 will end up with sim- 
ilar prospects, whereas the pseudo-scalar O2 somewhat 
worse prospects, because it leads to a longer \* lifetime, 
with decays typically happening outside of the detector. 
In that case, the standard monojet search will probably 
better probe the iDM model. 

Conclusions. In conclusion, we have examined iDM 
models whose larger splitting precludes their direct de- 
tection (at tree level - scattering at one loop level is 
nonetheless possible, but will occur with negligible rate 
for A - 1 TeV). We find that the LHC has the ability to 
search for such models for a wide variety of masses and 
splittings, though a novel signature involving a monojet 
plus extra hadrons which are reminiscent of hadronic tau 
decays, but appear deep in the detector. We hope these 
initial promising results inspire more detailed analyses by 
the experimental collaborations. 
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